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O Abstract 

^< ', The equivalence between f(R) gravity and scalar-tensor theories is invoked to study the null, 

• strong, weak and dominant energy conditions in Brans-Dicke theory. We consider the validity 

of the energy conditions in Brans-Dicke theory by invoking the energy conditions derived from 
a generic f(R) theory. The parameters involved are shown to be consistent with an accelerated 
O ', expanding universe. 
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CN ■ 1 Introduction 

Recent observations have revealed that the present state of the universe is undergoing an accelerated 
expansion [1]. There are, in general, a number of different approaches towards explaining this accel- 
eration. One such approach utilizes what is known as f(R) modification of gravity which, in effect 
is equivalent to Brans-Dicke (BD) type theories. The assumption of the existence of dark energy 
is another approach often used in this respect. In all the above theories, energy conditions impose 
stringent constraints whose validity should be studied in the light of their ability to explain the obser- 
vational data. From a theoretical viewpoint, energy conditions in their various forms, namely strong 
energy condition (SEC), weak energy condition (WEC), dominant energy condition (DEC), and null 
energy condition (NEC) have been used in different contexts to derive general results that would hold 
for a variety of situations [2]. For example, the Hawking-Penrose singularity theorems invoke the 
WEC and SEC [3], whereas the proof of the second law of black hole thermodynamics requires the 
NEC [4]. Another example comes from cosmology [5] where energy conditions are studied by using 
red shifts. 

The equivalence between BD theories and f{R) gravity is a subject that has been studied by 
various authors, as an example see [6]. The study of energy conditions may thus benefit from such 
equivalence in that knowing the energy conditions in one theory would point to the energy conditions 
in the other. For example, in the BD theory case V{4>) = Vocf 2 [7] where Vq cannot be found easily, 
one can use the equivalent f(R) theory to facilitate the calculation of Vq which, could then be used 
in the BD theory. 
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Even though the goal of this paper is to study energy conditions in modified theories of gravity 
and consequently in Brans-Dicke theory, much of the techniques will be borrowed from the analysis 
of energy conditions in Einstein's gravity. Therefore, we shall briefly review the derivation of energy 
conditions in general relativity. For a pedagogical review, see for example [3]. 

2 Energy conditions in general relativity 

Let v^ 1 be a tangent vector to a congruence of time like geodesies. For a hypersurface orthogonal 
congruence, Raychaudhuri's equation reads 

^ = -\Q 2 -°^ v -R^v\ (1) 

where 6 and a^ u are the expansion and sheer of two nearby tangent vectors, respectively. In Einstein's 
gravity, SEC and the ensuing singularities theorem follow from requiring that 

R^tfv" = 8ttG(T^ - X - 9iXV T)v»v v > 0, (2) 

for all time-like for which ^ < 0. A pair of nearby time-like geodesic vectors converge and 
will eventually intersect. The stress-energy tensor at each point p 6 M = K 4 obeys the inequality 
T^ u v^v u > for any time like vector v M G T p of an observer whose world-line at p has the unit 
tangent vector v and the local energy density appears to be T^ v . This assumption is thus equivalent 
to the energy density being non-negative as measured by any observer which, of course, is physically 
reasonable. 



3 Energy condition in f(R) gravity 

In this section we use the metric formalism in f(R) gravity and derive the strong, weak and dominant 
energy conditions for a general form of f{R). In doing so we will follow the formalism recently 
developed in [8]. We take the Freedman- Robertson- Walker (FRW) metric to study the cosmological 
implications of the models studied here. 
The action for f(R) gravity is [9] 

S = j^j d A x^f(R)+S m , (3) 

where we have set n = 8irG = 1. The field equations resulting from this action in the metric approach, 
assuming the connections are that of the Levi-Civita, are given by 

Gjxv = R^v ~ ~^9jxvR = T^ u + T^ v , (4) 

where T™ = jtj^ represents the energy-momentum tensor of ordinary matter considered as perfect 
fluid given by 

T 'lTu = (P + P)u^u v + pg^, (5) 
and T° v is the stress energy tensor of the gravitational fluid 
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l 9Af{R) - Rf'(R)) + (g^gpu - 9^g a/ 3)V a ^f'(R) 



» v f'(R) [2^ 

where a prime represents differentiation with respect to R. The field equation (4) now reads 



(6) 



f\R)R^ - \f(R)g^ - (V„V„ - g^a)f'(R) = T™. (7) 
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Contracting the above equation we obtain 

f'{R)R - 2f(R) + Saf(R) = T m . 



(8) 



Now, let us briefly review the energy conditions in f(R) gravity. We begin by defining an effective 
stress-energy tensor using equation (7) as follows 



peff _ 
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T™ + ~ {f(R) - Rf'(R)) 9»u + (V M V„ - g^a)f(R) 



with 



T cS = jy^r [T m + 2(f(R) - Rf'(R)) - 3Bf(R)] . 



(9) 



(10) 



Therefore, we can write R^ in the terms of an effective stress-energy tensor and its trace, that is 



Using the spatially flat (k = 0) FRW metric as 

ds 2 = -dt 2 + a(t) 2 dx 2 , 
the effective energy density and pressure are given by 



neflf 



„eff 



1 



1 



p--(f(R)-Rf'(R))-3HRf"(R) 



and 



P 



eflf 



f(R) r 2 



f(R) r 2 



p + -(f(R) - Rf'(R)) + (2RH + R)f'(R) + R 2 f"'(R) 



(11) 



(12) 



(13) 



(14) 



where R = dR/dt and H{t) = is the Hubble parameter. Now, using these equations, we can 
write the NEC and SEC, given by p cS + p cS > and p cS + 3p cS > respectively as 

p + p + (R - RH)f'(R) + R 2 f"'{R) > (15) 

and 



P + 3p + (/(i?) - Rf'(R)) + 3(R + RH)f"(R) + 3R 2 f"'(R) > 



(16) 



To compare our results here with that of general relativity for a given f(R), we use the FRW 
metric which, for WEC (p cS > 0) leads to 



p - -(f(R) - Rf'(R)) - 3HRf"(R) > 0. 



(17) 



For DEC (p etl - p etl > 0), we find 



p- p - (§RH + R)f"(R) - R 2 f"'{R) - (f(R) - Rf'(R)) > 0. 



(18) 



4 Energy conditions in Brans-Dicke theory 



Let us now investigate a non-minimally coupled self interacting scalar-tensor field theory such as the 
Brans-Dicke (BD) theory and find the various energy conditions for this type of modified gravity. In 
the context of BD theory [10] with a self interacting potential and a matter field, the action is given 
by 



xv-g 



<j>R - -V a <j>V a <p - V{<t>) 



(19) 



3 



where u is the usual BD parameter and we have chosen units such that 8irG = c = 1. The gravitational 
field equations can be derived from action (18) by varying the action with respect to the metric 



G 



- flV 



2(f> ' 



(20) 



where T™ v is the stress-energy tensor of the normal matter as expressed in equation (5). Variation of 
action (19) with respect to 4> gives 



n<f> 



T m 1 
+ 



dV{4>) 



2uj + 3 2^ + 3 
where the expression n<j> = g^ u V ^V u (j), using (12) is given by 



2V{4>) 



(21) 



(22) 



Using the equation of motion, we can write = k'{T™ + T$ u ) where k' = ^. The stress-energy 
tensor and its trace for the BD theory may now be calculated with the result 



and 



= t[(V</>) 2 - 2V a (j)V a (p} - 3D<j> - 2V{6). 
<P 

Comparison of these equations with equation (2) leads to a similar equation for the SEC, 



R^v u = k' 



(T™+T*)--g, u (T m + T*) 



> 0. 



(23) 
(24) 

(25) 



We can now write the relations for NEC and SEC analogously as (p + p^ + p + p^) > and (p + p^ + 
3{p+p^)) > respectively [3] so that we may first derive p^ and p^ for the spatially flat FRW metric 
as follows 

,,, . V(rh\ 

(26) 



2/ -?> H <i> + — 



and 



P* = fy 2 + (2H<j> + . 
Thus, the NEC and SEC for the BD theory are given by 

P + P+^r<p 2 + (4>- H<j>) >0, 



(27) 



(28) 



and 



2uj -o 

p + 3p+ —<j) 2 + 3(0 + H<j>) - V{4>) > 0. 



(29) 



Now, following and expanding on the GR approach to include WEC and DEC, as has been employed 
in f(R) gravity theories, we may obtain similar equations in the BD theory. Therefore, the WEC 
and DEC in the BD theory are respectively given by 



and 



P-P- (4> + 5H<j)) +V(<f>) > 0. 



(30) 
(31) 
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5 Equivalence of the energy conditions in f(R) gravity and Brans- 
Dicke theory 

Considering action (3) within the context of the metric formulation of f(R) gravity, one can introduce 
a new field \ and write a dynamically equivalent action [6] 

S= hj dAx ^9[f{x) + f'{x){R-x)] +«Sm(W)- (32) 

Variation with respect to x leads to equation % = R provided f"(R) / 0, which reproduces action 
(3). Redefining the field x by = f'(R) and setting 

V( ( f>) = M^-f(xm, (33) 

the action takes the form 

S = 7 LJ dAx ^9^ R - V ^ + S m(9^^). (34) 

Comparison with action (19) reveals that this is the action of a BD theory with lo = 0. Therefore, 
metric f(R) theories, as has been observed long ago, are fully equivalent to a class of BD theories 
with a vanishing kinetic term [6]. Now, taking x = R and <f) = f'(R) we have 

4> = Rf"(R) and cf> = Rf"(R) + R? f"'(R). (35) 

Substituting these relations into equations for NEC, SEC, WEC and DEC in BD theory, namely 
equations (28), (29), (30) and (31) with lo = , one can easily derive respectively NEC, SEC, WEC 
and DEC for the f(R) modification of gravity, that is equations (15), (16), (17) and (18). 

5.1 Examples 

To see how equation (17) can be used to put constraints on a given f(R) and equivalently on the BD 
potential, let us examine two examples. First, we consider f(R) as having a general power-law form, 
given by 

f(R) = oR n . (36) 

Let us now concentrate on the vacuum sector i.e. p = p = 0. Substituting in equation (17) we have 
the following condition for WEC 

a(n- 1)(1 -nA) > 0, (37) 

where A = _ g ° J± and the deceleration (q ), jerk (j ) and snap (s ) parameters for the present-day 

values are defined in [8]. I what follows, we examine two values of the exponent n, namely n = — 1 
and n = 2 which satisfy the inequality (37). 

Case I: n = — 1 

Takeing q ~ —0.81 and j ~ 2.16 given in [8], so that A ~ 0.29. Equation (37) for n = — 1 reduces 
to 

-2.58a > 0. 

This relation is satisfied for a < which in turn satisfies equation (37) with n = — 1. The potential in 
BD theory with lo = for V((f>) = Vo4> m , corresponding to f(R) = crR n can thus be obtained simply 
as 

_ (n-l)a 
(na) m 

where m = so that for n = — 1 the corresponding BD potential is V((f>) = —2\J—a^h. 
Case II: n = 2 

As a second case we consider n = 2, so that f(R) = aR 2 , leading to a WEC given by — 0.06<r > 
which requires a < 0. The corresponding BD potential is V(<p) = j^4> 2 with a being negative. One 
can therefore come to the conclusion that Vq must also have a negative value. 
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6 Redshift and energy conditions 



Let us now take the potential V{4>) = Vo(j) m in BD theory with u = and use the following power- low 
ansatze 

tV . „ ft 



\) i — j and a{t) = a y— j ■ (38) 

Inserting these relations into equation (21) in the vacuum sector, one then finds that 

(5 = (39) 
1 — m 

Substituting the spatially flat FRW metric (12) in the field equations (4) we get 

»(!)'=/■, (40) 

and 

■« ^«^ 2 „eff 



2-+ - =-p ett . (41) 
a \a J 

Now, from equations (38) and (40) we can write 

p eS = 3tf 2 (l + z)i (42) 

where we have used the relations = 1 + z, Ho = ^|^y and a = to-^o where z is the redshift of a 
luminous source [11]. From this equation one may conclude that in the very early universe a would 
have been small and therefore p eS large. Now, the WEC is given by 

3flg(l + z)« > 0. (43) 

This equation for the effective energy density clearly satisfies the WEC. If we follow the same method 
as in section 3 or in [8] and substitute the total energy density and pressure by p eS and p eS , we will 
find the same relation between the distance modulus and redshift parameter which is studied in [11] 
where the energy conditions have been used. Now let us write the NEC, DEC and SEC with respect 
to the redshift parameter as follows 

2-Ho,- ^ 



to 



-(l + z)->0, (44) 



(6a- 2) — (1 + z)i > 0, (45) 
to 

6(1 -a) — (l + z)« >0, (46) 
to 

respectively. From equation (44) one can see that it is not possible to extract more information from 
NEC. As far as the DEC and SEC are concerned however, they require a > | and a < 1 respectively. 
The first is in agreement with the observation that the universe is undergoing an accelerated expansion 
phase, a > 1. A glance at the last equation reveals that it is completely in contradiction with an 
accelerated expanding universe, but we know that the SEC ensures gravity to be always attractive. 
Violation, as discussed in [12], allows for the late time accelerated cosmic expansion as suggested by 
the combination of recent astronomical observations. 
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7 Energy condition in the Einstein frame 



What we have done so far in the pervious sections has been in the so-called Jordan frame. However, 
it would also be instructive to study these relations in the Einstein frame. As is well known, the 
usual procedure, going from one frame to the other, is to use a conformal transformation. A problem 
then arises in that whether the tensor representing the physical metric structure of space-time is the 
one belonging to the Jordan frame or to the Einstein frame. However, what we are concerned with 
in this work is the relation between the energy conditions in these frames and will not deal with the 
question posed above. Under the conformal transformation 



9^ e^V, (47) 
<f> = -In f'(R), (48) 



and taking 
the action (3) is rewritten as 

S E = J d 4 x^9 (R ~ \g^d a ^ - V(<P^j 



(49) 



Here 

R /(g) 
f(R) f'(R) 

As a result, one finds the field equations for the metric in the form 



v(4>) = - 4rnk ■ (50) 



G,xv = |v„0V„0 - ^g^g af3 V a (f>Vp4> - ^V(^ u . (51) 



Also, the equation of motion for <f> becomes 
For the FRW metric we find 

n</) = - e - <t '(4> + ^ 2 + 3H<})). (53) 
In order to solve these equations for the case f(R) = crR n we use the following power-low ansdtze 

a(t) oc P and <f> = alnt, (54) 

where a and (3 are arbitrary constants. The Hubble parameter, H, then reads 

H = i . (55) 



Now, using equations (48) and (50) we have 

V( , ) = ^l(^f i . (56) 

Matching the exponents of t in equation (51) we arrive at the following expression for a 

2(1 -n) 

a = -2^Y- (57) 



If we define Tfi v as 



= \v»<j>V v <j> - l~g^9 a ^ a <pV^ - \v{<f>)~g^, (58) 
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then we obtain the following relations for and 



+ -e<V(0) and = ^LZ e^V^a 2 . (59) 



4 2 "12 
Now let us write the WEC, SEC, DEC, and NEC 

3<j) 2 



2 +e^)>0, (60) 



^l±|^! + e ^)(l-3a 2) >0, (61) 



,,2 ,2 

2~ 

,2\i2 



fcfWL + eM«(l + o 2 )>0, 



2 

Sli^! + eV W (l-^)>0. 
Using our ansdtze (54), we find the following expression for the WEC 



(62) 
(63) 



2— n 

6(l-n) 2 /n-l\ / 1 \~ , t . 

As we can see in equation (64), the WEC always holds. However, in the case that either n < or 
a < 0, we have a constraint on n, namely that n must be even. For the other energy conditions at 
late times, a(t) — ► oo, we have 

2 2 — n 

71 \ i 71 — 1 \ / 1 \ n ~ 1 

" ' ' ^ ' ' - J > SEC and NEC, (65) 



2n — 3 / \ " 2 <7 / \na J 

Equation (65) tells us that there is no guarantee that it remains positive. As can be seen, when the 
SEC and NEC hold, the DEC does not. For these equations we also have a constraint on n, that is, 
when either n < or a < 0, then n must be even. In the case of n < another constraint appears 
which is important for both equations. There are ranges of n for which either SEC or NEC do not 
hold. Figure 1 shows this situation. Also, for the case n > 0, we see that the DEC does not hold at 
all, regardless of the sign of a. 



8 Conclusions 

In this work we have studied the energy conditions in the BD theory and compared the results with 
that of f(R) gravity, benefiting from the ease with which the parameters of interest can be derived 
in the latter and subsequently used in the former. This would help us to check the validity of the 
energy conditions in BD theory by invoking the energy conditions derived from a generic f{R) theory. 
The parameters involved were shown to be consistent and compatible with an accelerated expanding 
universe. 
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Figure 1: The behavior of SEC and NEC as a function of n for n < and a > 0, left and of SEC and NEC for n < 
and a < 0, right. 
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